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Abstract 

We present a full next-to-leading order (NLO) QCD corrections to diphoton pro- 
duction at the hadron colliders in both standard model and ADD model. The invariant 
mass and rapidity distributions of the diphotons are obtained using a semi-analytical two 
cut-off phase space slicing method which allows for a successful numerical implementa- 
tion of various kinematical cuts used in the experiments. The fragmentation photons 
are systematically removed using smooth-cone-isolation cuts on the photons. The NLO 
QCD corrections not only stabilise the perturbative predictions but also enhance the 
production cross section significantly. 



mc.kumar@saha.ac.in 
^prakash . mathews@saha. ac . in 
■^ravindraQmri.ernet.in 
^ anurag@mri . ernet . in 



1 Introduction 



The gauge hierarchy problem has been one of the main motivations to go beyond the 
standard model (SM). A novel idea that addresses this problem was put forward by 
Arkani-Hamed, Dimopoulos and Dvali (ADD) wherein they introduced extra spatial 
dimensions and allowed only gravity to propagate in the extra dimensions, keeping the 
SM fields confined to a 3-brane [1] . As the inverse square law behavior of gravity has so far 
been tested down to sub-millimeter length scales, the size of the extra dimensions, in this 
model, should be much smaller than sub-millimeter. The apparent weakness of gravity 
as compared to the other forces seen in nature, can now be accounted for through the 
volume of the extra dimensions. The relation between the fundamental scale Ms at which 
the new physics sets in (above which the extra dimensions are dynamically accessible) 
and the Planck scale Mp is given by 

Ml ~ Mf^''^ , (1) 

where d is the number of extra spatial dimensions and the size of the extra dimensions. 
Since R is of order of a milli-meter, the scale Mg can be as low as a few TeV, which 
circumvents the hierarchy problem. The propagation of a massless graviton in 4 -|- c? 
dimensions, after compactifying the extra dimensions on a rf- dimensional torus, manifests 
itself as an infinite tower of massive Kaluza-Klein (KK) modes on the 3-brane. Each 
KK mode couples with SM field through energy momentum tensor with a coupling 
proportional to k ~ 1/Mp. However, the effective coupling after summing over all the 
KK modes is enhanced significantly due to large multiplicity of KK modes. In any typical 
scattering process at colliders, the gravity can enter through their KK propagator as well 
as through the real emission of KK states. These KK states are large in number. Hence 
the suppression resulting from coupling k is compensated by the large multiplicity factor 
resulting either from the sum of KK propagator T>{Q'^) or from the phase space of large 
number of real KK states . For example, if the KK states enter through a propagator. 
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we find that any typical amplitude will be proportional to 

n " 

^ H^ + 2/(A/Q)l, (2) 



where A = is the explicit cut-off on the KK sum and the function I can be found 
in [2]. Thus, for Mg ~ O(TeV), the gravity effects can become significant and hence the 
collider phenomenology associated with this model is very interesting [2]. To exemplify, 
the virtual effects of the KK modes could lead to the enhancement of the cross sections 
of pair productions in the processes like Drell-Yan, diphoton and dijet while the real 
emissions could lead to large missing signals giving some new observable like mono- 
jet, mono-photon in an experiment. Owing to a very high centre-of-mass energy of 

= 14 TeV and a large gluon flux at the large hadron collider (LHC), rich collider 
signals resulting from this model have been reported in the literature [2-7]. However, 
these results are based on leading order (LO) calculations. At the hadron coUiders hke 
LHC, the QCD effects are often considerably large and hence the quantum corrections 
can influence the predictions significantly. In the ADD model, QCD effects [8] have 
been shown to increase the di-lepton productions and also to stabilise the perturbative 
predictions. Hence in this paper we study the impact of the QCD corrections for the 
diphoton signal in the ADD model. 

In QCD, the infra-red safe observable exhibit a feature called factorisation, according 
to which collinear singularities can be factored out from the partonic cross sections in a 
process independent way and then they are either absorbed into the bare parton distri- 
bution functions (PDF) if they originate from initial state partons or into fragmentation 
functions if they are from final state partons. This procedure introduces a scale called 
factorization scale fip, which is arbitrary. In addition, ultra- violate renormalisation in- 
troduces renormalisation scale Hr which is again arbitrary. The truncated perturbative 
expansion leaves our theoretical predictions /ip and /ir dependent, these scale depen- 
dence will go down as we include more and more terms in the expansion. In addition. 



the fitted PDFs are usually not fully constrained due to insufficient experimental data. 
Hence, predictions beyond LO are often more reliable than LO ones. 

Diphoton production process is an important probe for the Higgs boson search at 
the LHC NLO QCD corrections to this process in the SM are available in the literature 
[9-12] and hence the diphoton signal has been a useful tool for precision studies. This 
process has also been used to search for the physics beyond the standard model, such 
as extra dimensional models, super symmetry and the unparticle physics. Di-photon 
production [5] at Tevatron has set stringent constraints on the parameters of the ADD 
model [13]. It will also play an important role at LHC. The D0collaboration [13] assumed 
a K-factor for their analysis but a full NLO QCD calculation for the ADD model does 
not exist for the diphoton production. In this paper, we have systematically computed 
all the QCD effects to NLO in perturbation theory to various important observable in 
diphoton production that are sensitive to the ADD model. Quantitative estimates of 
QCD corrections to these observable are presented and our predictions are expected to 
be less sensitive to the factorisation scale. 

2 The Diphoton Production 

In the SM, at leading order (LO), diphoton production proceeds via quark anti-quark 
annihilation subprocess g + g^7 + 7u. In the ADD model, the SM fields couple to KK 
modes through the energy-momentum tensor of the SM fields with a strength denoted by 
K. Hence, diphotons are produced in (i) quark antiquark annihilation (5'+^ ^ 7 + 7) and 
(ii) gluon fusion process {g + g — >■ 7 + 7) via the exchange of KK modes. A comprehensive 
phenomenology taking into account all the above LO processes has been done in [5]. It 
was observed that unitarity restricts the maximum value of the invariant mass Q of the 
diphotons. Following [5], we restrict the invariant mass Q to Q < 0.9 Mg. 

At NLO, the SM as well as ADD leading order quark antiquark annihilation processes 



^The gluon-giuon fusion process through quark loop, though of order a^, is comparable to the LO 
for studies of photon pairs having small invariant masses, M^-y. As it falls off rapidly as M^^ increases, 
it no longer enjoys the status of LO process for our study on the production of large invariant mass 
photon pairs in the context of ADD model and is truly a NNLO contribution. 
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get 0{as) QCD radiative corrections through virtual gluons ing + g— ^7 + 7 + one loop 
and real gluon emissions in q+q 7+7+5' processes. To this order, q{q)+g Q'(^)+7+7 
process also shows up in both SM and ADD. The LO gluon fusion process in the ADD 
model gets NLO QCD corrections to order ag through fl' + g'— >7 + 7 + one loop and 
g + g^^ + ^ + g processes. Since KK modes appear at the propagator level, the LO 
SM (ADD) processes interfere with the corresponding NLO ADD (SM) processes giving 
order ag NLO QCD corrections. We have incorporated all these NLO QCD corrections 
in this article for the study that follows. 

The NLO partonic cross sections are often ill-defined due to soft and coUinear singu- 
larities that result from the presence of zero momentum gluons and mass-less partons. In 
addition to these singularities, we encounter coUinear (QED) singularities that originate 
when the photon in the final state becomes coUinear to the quark or the anti-quark emit- 
ting it. These (QED) singularities go away if we also include the diphoton production 
channels resulting from the fragmentation of partons. This involves introduction of non- 
perturbative fragmentation functions. These functions are poorly constrained. Hence, in 
our study we do not include fragmentation photons but consider only direct photons. Al- 
ternatively, we can suppress QED coUinear singularities using the smooth-cone-isolation 
prescription proposed by Frixione [14]. In the rapidity-azimuthal angle {y,(p) plane the 
amount of transverse hadronic energy Et in any cone of radius r = -y/ (Ay)^ + (A0)^ 
with r < To centered around the photon must satisfy 



The above prescription safely removes all the photons from the fragmentation processes 
without disturbing soft and coUinear partons. 

An analytical computation incorporating smooth-cone-isolation and other kinemat- 
ical constraints at NLO level is hard to achieve. Hence, we resort to a semi-analytical 
approach called two cutoff phase space slicing method [15]. In this method, two small 
slicing parameters Sg and 6c are introduced to isolate the cross sections that are sensitive 




(3) 
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to soft and collinear singularities. The remaining part of the cross section denoted by 
da^'^"'{6s, Sc) is soft and coUinear free. The soft divergences come from virtual as well as 
real gluons when their momenta become zero. On the other hand the collinear singular- 
ities arise due to mass less nature of the partons. We compute these soft and collinear 
sensitive cross sections (they are singular in 4 dimensions) analytically in 4 + e dimen- 
sions which regulate these singularities. The soft singularities cancel between virtual and 
real gluons when their contributions are added appropriately. The remaining collinear 
singular terms which appear as poles in e are systematically removed by collinear counter 
terms in MS factorization scheme. This is usually done at an arbitrary scale fip. Hence 
we will end up with a finite cross section coming from (a) soft and collinear sensitive 
regions denoted by d(7''^''^*"'(5s, 5c) {sc denotes soft and collinear) and (b) (i(T-^*"'(5s, 5c) 
part of the cross section. Their sum, ie. (a) + (6), is expected to be free of choice of the 
slicing parameters. This is an essential prerequisites for the implementation of the phase 
space slicing method. 

3 Numerical Results 

In this section, we present our results for invariant mass (Q) and rapidity (Y) distri- 
butions of the photon pair at LHC We have employed the kinematical cuts given by 
ATLAS collaboration [16]: the transverse momentum > 40 GeV for the harder pho- 
tons, p'^ > 25 GeV for the softer photon, and the rapidity \y^\ < 2.5 for each photon. 

In addition, the photons are isolated from hadronic activity according to Eq. ([3]), with 
n = 2,ro = OA,E^° = 15 GeV. The minimum separation between the two photons is 
taken to be r^^ = 0.4. For the LO, we have used CTEQ6L PDFs and CTEQ6M for NLO 
studies [17], with the corresponding value of a<j(M^) = 0.118 and 5 light quark flavours. 
The factorisation and renormalisation scales are taken to be Q, the invariant mass of the 
diphoton pair. The electromagnetic coupling constant is chosen to be a = 1/128. For 
our numerical analysis we have chosen the ADD parameters, Mg = 2 TeV and A = M<. 
for the number of extra spatial dimensions d = 3. This choice of = 2 is consistent 
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Figure 1: Stability of the order as contribution to the total (SM+ADD) cross section 
against the variation of the phase space slicing parameters 6s (left) and Sc (right) in the 
invariant mass distribution of the di-photon system with Ms = 2 TeV and c? = 3 at 
Q = 700 GeV. 



with the limits from [13]. 

We have first checked our numerical code by studying the dependence of observable 
on the slicing parameters, 6s and 6c- In the left (right) panel of Fig. [T] we have plotted the 
order contribution to the invariant mass distribution of diphotons in SM and ADD 
against the slicing parameter 6s {6c) in the range between 5 x 10""^ and 10~^. For the 
6s variation (left panel) we have fixed 6c = 10~^ and for the 6c variation (right panel) 
we have fixed 6s = 10^^. These plots show that our numerical results are least sensitive 
to the slicing parameters for a wide range. The percentage of uncertainty that results 
from the choice of slicing parameters is found to be around 6.7%. This study confirms 
the reliability of our code for further predictions. For our numerical predictions, we have 
chosen 6c = 10~^ and 6s = 10^^. Other important check on our code comes from a 
detailed comparison of our SM results against those given in the literature [9-12]. In 
particular, we find that our SM results are in very good agreement with those given 
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Figure 2: Various subprocess contribution to tlie invariant mass distribution of the dipho- 
ton production with = 2 TeV and d = 3. The SM gg subprocess (lower solid line) is 
at 0{a'^) while all other subprocess are at order 0{as)- 

in [11] with their choice of parameters. 

In Fig. [21 we have presented various subprocess contributions to the invariant mass 
distribution of the diphoton system for the range 400 < Q < 1100 GeV where gravity 
(through KK modes) contribution dominates over the SM. Both qq and gg initiated sub- 
processes in ADD give large positive contributions while the qg initiated subprocess gives 
a negative contribution. The interference of the SM with ADD (SM*ADD) from both 
qq and gg subprocesses gives almost Q independent contribution, while the contribution 
from the qg subprocess falls steeply at higher values of Q. Owing to the large gluon flux 
at the LHC, the gg initiated subprocesses in ADD give the dominant contribution over 
the rest, thus making the observable effects of ADD model clearly visible in the large 
Q region. We have also plotted the SM gluon-gluon fusion sub-process through quark 
loop contribution separately in the Fig. [21 It is clear from the plot that its contribution 
is negligible compared to SM quark anti-quark initiated processes and hence belongs to 
NNLO contributions. Hence, we have not included this in our study. 
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Figure 3: Invariant mass (left) and rapidity (right) distributions of the diphoton pro- 
duction at the LHC with Mg = 2 TeV and d = 3. For rapidity distribution, we have 
integrated over Q in the range 600 < Q < 1100 GeV. 

In Fig. m we have presented the invariant mass (left panel) and rapidity (right panel) 
distributions of the diphoton productions in both SM and ADD model. We have plotted 
LO and NLO contributions separately to demonstrate the impact of QCD corrections. 
It is clear from the plots that the QCD corrections to both invariant mass and rapidity 
distributions in SM as well as in ADD model are large for the entire range of Q considered. 
In the left panel we find that the contribution from ADD dominates over that of SM 
starting around Q = 500 GeV. The exact value where this happens depends crucially 
on the parameters of ADD model. For the rapidity distribution (right panel), we have 
considered \Y\ < 2.0 and integrated over Q in the range 600 < Q < 1100 GeV where 
the KK effects are dominant. The cross section is found to be maximum at the central 
rapidity region both in SM and in ADD model, the later differing by more than an order 
of magnitude. 

The cross sections do depend on the isolation criterion. The at the partonic level 
need not be the same as that of the hadrons at the detector level, which gives rise to 
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Figure 4: Dependence of the invariant mass distribution of the diphoton system on the 
parameters (left) and n (right), of the Frixione's isolation algorithm, with Mg = 2 
TeV and d = 3. For the variation of E?^° (n) we have kept n {E^°) fixed. 

the dependency of the cross sections on E^°. In the smooth cone isolation prescription 
discussed above, large logarithms of often spoil the reliability of fixed order compu- 
tation. We can study the effect the these logarithms by varying the function that appear 
in the isolation criterion. We present in Fig. HI the dependency of our cross sections 
on the choices of E^° (varied between 5 GeV and 30 GeV), and n, (varied between 1 
and 2). We find that the dependency is unnoticeable making our predictions reliable for 
experimental study. 

Finally we consider the invariant mass distribution at the Tevatron for both the SM 
and ADD model to NLO QCD. We have used Ms value which is consistent with the 
experimental bounds [13] for the di-electromagnetic signal which is the combined e~^e~ 
and 77 final state. In this analysis we are hence interested only in gauging the impact of 
the QCD corrections to these studies. In Fig. [5] we plot the invariant mass distribution 
of the diphoton system in the range 100 < Q < 1000 GeV at the Tevatron (^/S = 1.96 
GeV) for both the SM and including the ADD contribution at LO and NLO in QCD. 
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Figure 5: Invariant mass distribution of the diphoton at the Tevatron for Mg = 2 TeV 
and d = 4 and in right panel the various contributing subprocess. The SM gg subprocess 
is the lower solid line which is at while all other sub process is plotted at 0{as)- 

We have used the following kinematical cuts: (a) transverse momentum > 15 (14) 
GeV for the harder (softer) photons, (b) rapidity \y^\ < 1.1 for each photon, and (c) 
vq = 0.4 and r^^ = 0.4. In addition for the smooth-cone-isolation we use = 2 GeV 
and n = 2. The contributions of the various subprocess is shown in the right panel, for 
the range 400 < Q < 1000 GeV. We have used the number of extra spacial dimensions 
d = 4 and Mg = 2 TeV. The impact of QCD corrections at the Tevatron is much mild 
compared to the LHC where the gluonic flux is overwhelming. 

4 Conclusions 

In this article, we have systematically computed NLO QCD corrections to the diphoton 
production process at the hadron colliders in SM as well as in ADD model. We use a semi- 
analytical two cut-off phase space slicing method to compute invariant mass as well as 
rapidity distributions of the diphotons system. We have applied the kinematical cuts used 
by the ATLAS detector collaboration for our study. A smooth-cone-isolation prescription 
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on the diphotons has been used to reject poorly known fragmentation photons. Our 
method takes care of all the soft and coUinear singularities that appear at NLO level in 
QCD. We have explicitly shown that our NLO results are least sensitive to the slicing 
parameters Ss and Sc- Our SM results are in good agreement with those given in the 
literature. Predictions for invariant mass distribution of diphotons in ADD model with 
Ms — 2 TeV are found to be large compared to those in SM for invariant mass Q > 600 
GeV. This is due to large gluon flux at the LHC which enhances the ghion initiated 
production channels over the rest. In addition, the QCD corrections are significantly 
large both in the SM and in the ADD over the entire range of Q considered. For the 
rapidity distribution, we have integrated Q in the region 600 < Q < 1100 GeV where the 
gravity (through KK modes) contributes significantly. We find that the QCD corrections 
are important throughout the region |y| < 2.0. In addition, our results arc expected to 
be less sensitive to the uncertainties coming from the choice of factorisation scale. 

In summary, we have accomplished an important task of computing all the partonic 
contributions at NLO level in QCD to diphoton production at hadron coUiders both in 
SM and ADD model. These QCD corrections for the ADD model and its interference 
with the SM are being presented for the first time, while to this order the SM results 
already exist in the literature. The NLO QCD effects are found to be large and they 
are expected to reduce theoretical uncertainties, thus providing an excellent opportunity 
to put stringent bounds on the parameters of the ADD model when the experimental 
results are available. Quantitative impact of the NLO QCD corrections to both the ADD 
and RS model would be addressed in a future publication [18]. 
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